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SUMMARY 


This  study  focused  on  characterizing  matrix  cracking  and  delamination 
behavior  in  multidirectional  laminates.  Static  tension  and  tension-tension 
fatigue  tests  were  conducted  on  two  different  layups.  Damage  onset, 
accumulation,  and  residual  properties  were  measured.  Matrix  cracking  was  shown 
to  have  a  considerable  influence  on  residual  stiffness  of  glass  epoxy  laminates, 
and  could  be  predicted  reasonably  well  for  cracks  in  90°  plies  using  a  simple 
shear  lag  analysis.  A  fracture  mechanics  analysis  for  the  strain  energy  release 
rate  associated  with  90°  ply  matrix  crack  formation  was  developed  and  was  shown 
to  correlate  the  onset  of  90°  ply  cracks  in  different  laminates.  The  linear 
degradation  of  laminate  modulus  with  delamination  area,  previously  observed  for 
graphite  epoxy  laminates,  wa3  predicted  for  glass  epoxy  laminates  using  a  simple 
rule  of  mixtures  analysis.  The  strain  energy  release  rate  associated  with  edge 
delamination  formation  under  static  and  cyclic  loading  was  difficult  to  analyze 
because  of  the  presence  of  several  contemporary  damage  phenomena. 
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Inverse  length  parameter 

Lamina  Poisson  ratio 


INTRODUCTION 


The  application  of  composite  materials  in  modern  helicopters  is  being 
extended  to  more  complex  and  safety-critical  components  which  require  a  deeper 
understanding  of  the  material  behavior  to  guarantee  reliability.  In  the  future, 
civil  and  military  authorities  may  require  a  failsafe  or  damage-tolerant  design 
of  composite  helicopters.  This  need  is  recognized  in  principle  by  the  helicopter 
manufacturers,  but  it  is  clear  that  material  properties,  test  methods,  and 
application  of  composites  to  real  components  must  be  further  investigated  in 
order  to  achieve  reliable  structures.  Because  glass/epoxy  composites  are  already 
being  used  extensively  in  helicopter  rotor  systems,  an  investigation  of  the 
fatigue  behavior  of  glass/epoxy  laminates  was  undertaken  in  the  hopes  of  laying 
a  foundation  for  future  damage  assessment  of  helicopter  structural  components. 

Composite  laminates  commonly  employed  in  helicopters,  as  in  many  other 
applications,  usually  show  evidence  of  damage  before  failure.  The  usually 
adopted  "first  ply  failure"  criterion  [1j  can  very  roughly  predict  transverse 
ply  matrix  cracking  but  cannot  predict  delamination.  Both  of  these  two  damage 
modes  may  occur  early  in  the  life  of  the  component  but  may  have  only  a  small 
effect  on  overall  stiffness  and  strength  properties. 

The  scope  of  this  work  is  to  study  these  two  phenomena  from  a  fracture 
mechanics  point  of  view,  offering,  as  far  as  possible,  some  simple  but  reliable 
criteria  to  predict  onset  and  accumulation  of  matrix  cracking  and  delamination. 
Extensive  work  has  already  been  done  with  good  results  using  fracture  mechanics 
to  model  both  delamination  [2-b]  and  matrix  cracking  [5-13J,  but  only  in  a  few 
cases  are  these  results  readily  applicable  to  design. 


In  this  work  a  new  approach  is  developed  for  characterizing  matrix  cracking 

onset  and  growth.  The  critical  strain  energy  release  rate  for  formation  and 

accumulation  of  matrix  cracking  in  the  90°  plies  of  [0/90]  and  [±45/0/90] 
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glass/epoxy  composites  were  calculated  and  plotted  for  a  range  of  crack 
densities  in  the  transverse  plies. 

A  previously  developed  closed  form  model  is  available  to  predict  edge 
delamination  behavior  in  composite  laminates  [2].  In  the  present  study  this 
model  was  applied  to  glass/epoxy  quasi-isotropic  laminates,  where  the  presence 
of  pre-existing  matrix  cracking  was  observed.  The  theoretical  model  was  then 
compared  with  experimental  results.  Critical  strain  energy  release  rates 
associated  with  edge  delamination  were  measured  for  the  material  studied. 

EXPERIMENTS 

Two  E  glass/epoxy  laminates  were  tested  in  quasi-static  tension  and  in 
tension-tension  fatigue.  The  material  was  T1 90/X751  /50 ^ 1  ^  with  an  average  ply 
thickness  of  .225  mm,  50%  fiber  volume,  125  °C  curing  temperature.  The  measured 
tensile  failure  strain  in  the  fiber  direction  was  1.76?  and  the  following 
elastic  properties  were  measured  from  unidirectional  specimens: 

E,  =  43.7  GPa  E2  =  1 6 . 0  GPa 

G12  =  6.3  GPa  Vj  2  =  .284 

^Produced  by  Cyanamid  Fothergill  Ltd.,  U.K. 


The  layups  tested  were  L0/90]  and  [±45/0/90]  .  The  first  sequence  was 
chosen  to  isolate  the  90°  ply  matrix  cracking  from  other  damage  phenomena;  the 
second  sequence  was  chosen  to  highlight  the  edge  delamination  at  the  0/90 
interface,  because  of  the  high  peeling  stresses  which  are  generated  [2].  Tests 
were  run  on  an  Instron  servo-hydraulic  machine.  Axial  strain  was  derived  from 
displacements  measured  by  a  pair  of  LVDT's  over  an  82  mm  gage  length.  A  high 
intensity  light  was  shone  through  the  specimen  to  observe  and  photograph  the 
damage.  Matrix  crack  density  and  delamination  were  measured  from  high  contrast 
photographic  prints  like  the  one  shown  in  fig.  1.  For  each  laminate  five  quasi¬ 
static  tension  tests  were  run  with  a  ram  speed  of  0.5  mm/min.  At  specified  loads 
the  specimens  were  unloaded  to  record  damage  and  residual  stiffness  (tangent 
modulus).  Tension-tension  load-controlled  fatigue  tests  were  conducted  with  a 
repeated  sinusoidal  load  of  6  Hz  frequency  and  0.1  load  ratio;  20  fatigue  tests 
for  the  [0/90]  lay-up  and  25  tests  for  the  [+45/0/90]  lay-up  were  performed, 
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covering  the  whole  range  of  possible  fatigue  lives  from  1  to  106  cycles.  The 

cyclic  loading  was  interrupted  at  specified  intervals,  to  monitor  damage  and 

measure  static  stiffness. 

DAMAGE  OBSERVATIONS 


Figure  2  shows  schematically  the  damage  types  observed  on  quas i- i sotrop ic 
specimens,  namely  90°  ply  matrix  cracking,  edge  delamination  at  the  0/90 
interface  and  local  delaminations  originating  from  the  angle  ply  cracks. 
Depending  on  the  type  of  test  (  static  or  fatigue  )  and  on  the  load  levels, 
damage  accumulation  showed  some  qualitative  and  quantitative  differences. 

From  previous  works  it  is  known  that  matrix  cracks  do  not  grow  indefinitely 
in  number,  instead  they  reach  a  "saturation  value",  which  seems  to  be  a  function 


of  the  laminate  characteristics  [8j.  in  static  tests  failure  occurred  when  the 


number  of  90°  ply  matrix  cracks  was  still  low  and  angle-ply  cracks  had  not  yet 
developed.  The  edge  delamination  at  failure  was  extended  over  20%  of  the 
interfacial  area  and  the  stiffness  loss  was  typically  about  15%  of  the  initial 
stiffness. 

Fatigue  damage  accumulation  is  shown  in  fig. 3.  This  figure  shows  the 
results  of  a  single  test  run  at  160  MPa  maximum  sinusoidal  load  (R-0.1),  which 
is  about  60%  of  the  laminate  static  failure  stress,  and  is  typical  for  damage 
developed  during  fatigue  tests,  run  on  quasi  isotropic  glass  epoxy  laminates  . 

The  build  up  of  the  laminate  damage  is  plotted  against  fatigue  cycles. 
Figure  3(a)  shows  the  matrix  cracking  formation  in  90°,  +45°  (external)  and  -45° 
(internal)  plies.  For  each  of  these  plies,  matrix  cracks  reach  a  saturation 
level  at  different  number  of  cycles.  The  cracks  in  the  90°  plies  accumulate 
faster  than  matrix  cracks  in  the  others,  so  that  all  other  damage  may  be 
considered  to  happen  after  the  90°  plies  damage  reached  saturation.  Fig. 3(b) 
displays  the  edge  delamination  and  local  delaminations  growth.  Comparing 
fig. 3(a)  and  3(b),  the  delaminations  and  cracking  of  +45  and  -45  plies  appear  to 
be  almost  contemporary,  so  that  their  effects  must  be  superimposed.  The 
cumulative  effect  of  the  damage  on  the  specimen  longitudinal  modulus  E/E0  is 
shown  in  fig  3(c) . 

The  same  type  of  tests  were  performed  on  [0/90J  laminates  in  order  to 
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isolate  the  behavior  of  90°  ply  matrix  cracking.  The  L 0 / 9 0 j s  laminate  showed 
only  90°  matrix  cracking  in  static  tests.  Also  in  fatigue,  the  damage  was  mainly 
90°  matrix  cracking,  but  the  onset  of  matrix  cracking  in  0°  plies,  i.e. 
longitudinal  splits,  was  observed  before  final  failure.  Total  stiffness  loss 
just  before  failure  was  15%  in  static  tests  and  more  than  20%  in  fatigue  tests. 
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ANALYSIS  OF  MATRIX  CRACKING 


Matrix  cracking  occurs  when  tensile  stresses  normal  to  the  fibers  are 
present.  On  glass  reinforced  composite  materials  this  kind  of  damage  is 
responsible  for  a  considerable  fraction  of  the  total  stiffness  loss  both  in 
static  and  in  fatigue  tests.  To  predict  the  matrix  cracking  effect  on  different 
laminates  and  loadings,  available  models  have  been  studied  and  compared  to 
experimental  results,  and  a  fracture  mechanics  approach  has  been  applied  to 
characterize  matrix  crack  onset  and  growth. 

Stiffness  loss  in  a  [0  /90  J  laminate 

m  ns 

A  closed  form  expression  for  stiffness  I033  in  a  [0  /90  J  laminate  as  a 

m  n  s 

function  of  crack  density  D  (fig.  *4)  was  derived  previously  from  a  shear  lag 
analysis  [ 1 3 J  a3 
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.  +  nE,  tanh  (A/2D) 
1  mEt  A/2D 


(1) 


where 

3  G, ,  E„  (n+m) 

and  Ex,  E0  are  the  moduli  of  the  cracked  and  the  uncracked  laminate 
respectively.  The  Elt  E2  and  G12  are  the  lamina  longitudinal,  transverse  and 
shear  moduli  respectively,  and  h  is  the  ply  thickness. 
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To  evaluate  the  degree  of  approximation  introduced  by  eq.(1)  in  the 
laminates  presently  studied,  this  equation  has  been  compared  with  results  from  a 
NASTRAN  finite  element  model  (see  Appendix  A)  and  with  experimental  data.  Fig. 5 
compares  equation  (1)  (continuous  line)  and  Finite  Element  Model  (diamonds  and 
dashed  line)  with  data  obtained  in  static  and  fatigue  tests  on  [0/90] 

3 

laminates.  The  small  difference  observed  between  the  finite  element  results  and 
eq.(1)  will  produce  consequences  that  will  be  discussed  later.  The  comparison 
with  the  experimental  data  also  confirmed  the  validity  of  equation  (1)  for 
elevated  crack  densities. 


Stiffness  Loss  in  Other  Laminates 


Although  equation  (1)  is  only  valid  for  [CW90n]  stacking  sequences,  its 
use  can  be  extended  to  model  stiffness  loss  due  to  matrix  cracking  in  the  90° 
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plies,  surrounded  by  0°  laminae,  in  other  laminates  of  arbitrary  orientation. 
Previously,  a  technique  was  developed  for  estimating  the  reduction  in  E^  for  a 
cracked  90°  ply  as  a  function  of  the  crack  density  [6].  This  technique  is 
depicted  in  fig. 6.  First,  the  reduction  in  laminate  modulus,  E^/Eq,  of  a  [0/90]s 
laminate  is  determined  as  a  function  of  the  crack  density  in  the  90°  ply  (fig. 
6(a)).  In  the  current  study,  this  was  determined  using  eq.1  (fig. 5).  Next,  the 
linear  relationship  between  the  reduced  value  of  E 2  in  the  cracked  90°  ply  and 
the  C 0/90 J s  laminate  modulus  is  calculated  from  laminated  plate  theory  (fig. 
6(b)).  The  reduced  laminate  modulus,  E^/E^,  calculated  when  the  effective 
transverse  lamina  modulus  in  the  90°  ply,  E^/E^,  is  set  to  zero  is  designated 
E^/Eq,  and  is  the  value  obtained  from  the  ply  discount  method  where  the  cracked 
90°  ply  is  assumed  to  carry  no  load.  Finally,  the  equation  of  the  straight  line, 
relating  E^/E^  to  ex/eq  (fig<  6(b)),  is  combined  with  the  equation  describing 
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Ex/Eq  as  a  function  of  crack  density  (fig.  6(a))  to  generate  the  equation  for  a 
master  curve  relating  the  effective  transverse  modulus  of  a  cracked  90°  ply  in 
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an  arbitrary  laminate,  E£/E2>  to  the  crack  density  in  that  90°  ply  (fig  6(c)). 
This  master  curve  may  then  be  used  to  determine  the  appropriate  reduction  in 
for  a  cracked  90°  ply  in  another  laminate,  such  as  the  L46/-45/0/90J3  laminate 
tested  in  the  current  study. 

In  ref.  6  it  was  assumed  that  the  longitudinal  modulus  of  a  laminate  with 

matrix  cracks  is  a  linear  function  of  the  "  effective  damaged  lamina  " 

transverse  modulus  E'  normalized  to  the  undamaged  lamina  transverse  modulus  E2 

2 

(fig.  6(b)).  The  equation  of  this  straight  line  yields: 


E2 ' ( D , w)  E  (D ,w)  -  E 

-  =  JL - -  (2) 

E2  E0  -  E#  ' 

In  reference  6,  the  finite  element  method  was  used  to  determine  the  L0/90J 

3 

laminate  stiffness  loss  as  a  function  of  crack  density,  Ex(D,w),  whereas  in  the 
present  study  it  was  possible  to  start  from  a  closed  form  expression  (eq.1). 


& 


Equation  (a)  may  be  rearranged  to  yield 


E..  -  E*  ♦  (E0  -  EJ 


A  precise  expression  of  eq.  3  may  be  obtained  through  laminated  plate 


theory  (App.  B) ,  but  the  result  shows  no  significant  difference  from  eq.  3. 


Eq.  2  gives  E'/E2  as  a  function  of  crack  density  and  ply  thickness,  w=2nh 


(fig.  4),  for  a  laminate  stacking  sequence  where  the  relationship  between  Ex  and 


D*w  is  known.  In  the  present  case,  a  [_0/90Js  glass/epoxy  laminate  was  modeled 


using  equation  (1).  This  approach  is  valid  only  if  the  90°  ply  is  constrained  by 

adjacent  0°  plies.  For  other  cases,  where  more  complex  shear  stresses  are 

involved,  the  relationship  between  E^  and  crack  density  may  be  obtained  from 

finite  element  analysis  and  substituted  into  eq.2  to  determine  the  functional 

relationship  between  E],/E2  and  cracl<  density. 

To  check  the  validity  of  representing  the  effect  of  transverse  ply  cracking 

on  laminate  stiffness  in  terms  of  a  calculated  effective  modulus  E'2(D,w),  the 

longitudinal  modulus  of  a  [±45/0/90]  laminate  was  calculated  and  plotted  in 

s 

fig. 7  (solid  line)  for  several  90°  ply  crack  densities  using  E£(D,w)  for  the 
cracked  90°  plies  obtained  from  eq.2.  The  quasi-isotropic  laminate  containing 
various  crack  densities  was  also  modeled  using  the  NASTRAN  finite  element  code 
(App.  A).  The  laminate  stiffnesses  calculated  by  FEM  are  shown  in  fig. 7  as 
diamond  symbols  connected  by  a  dashed  line.  The  agreement  between  the  calculated 
stiffness  loss,  from  eq.2  and  FEM,  and  the  experimental  data  (triangles  in  fig. 


Strain  Energy  Release  Rate 


A  fracture  mechanics  analysis  was  applied  to  develop  a  tool  capable  of 
predicting  onset  and  growth  of  matrix  cracking.  For  an  elastic  body  containing  a 
flaw  of  area  A  that  extends  under  a  constant  nominal  strain,  e,  the  strain 
energy  release  rate,  G,  may  be  written  as 


e2  dE_ 
2  dA 


where  V  is  the  volume  of  the  body  and  dE/dA  is  the  rate  of  stiffness  change  as 
the  flaw  extends  [2]. 

The  fundamental  assumption  of  this  work  is  that  all  of  the  N  matrix  cracks 
in  the  volume  V  may  be  treated  as  a  single  equivalent  flaw  of  area  A,  that  is 


a  =  yN  a 
A  4-1  Ai 


and  that  all  the  matrix  cracks  have  the  same  length,  w,  equal  to  the  thickness 
of  the  cracked  laminae  (fig.  4).  Therefore,  for  a  laminate  of  width,  2b,  length, 
l,  and  total  thickness, t, 


A  =  N  2b  w,  and  V  =  2b  S,  t 


Because  N  =  i  D,  A  can  be  espressed  in  terms  of  crack  density  D: 


A  -  2b  1  w  D,  and  dA  »  2b  i  w  dD 


ms&rn. 


With  these  assumptions  G  may  be  written  as 


I 


i 


t  e2  dE 
w  2  dD 


Differentiating  eq.  1  and  substituting  into  eq.  7  yields  a  closed  form 

expression  of  the  strain  energy  release  rate  for  matrix  cracking  in  a  [0  /90  ] 

m  n  s 

laminate: 


tanh( X/2D ) 
X/2D 


cosh2 ( A/2D) 


(n+m)  E,  Ef 
m  E,  D 


r  +  nE,  tanh(A/2D) 

L  mEj  X/2D  J 


The  closed  form  expression  is  particularly  suitable  for  comparing  different 

lay-ups  and  materials.  Fig. 8  shows  the  calculated  effect  on  strain  energy 

release  rates  of  different  lay-ups  and  reinforcing  fibers.  In  the  case  of  Glass 

fibers  (solid  and  dashed  curves),  the  effects  of  doubling  either  the  90°  or  the 

0°  ply  thickness  are  shown.  Doubling  the  90°  ply  thickness  results  in  a  doubling 

of  the  strain  energy  release  rate,  as  it  is  evident  comparing  the  curves  for  the 

[0/90]  and  [0/90.]  laminates.  Doubling  the  number  of  0°  plies,  however,  has 
s  2  s 

very  little  effect  on  G  as  seen  when  comparing  the  curves  for  the  [0/90.]  and 

t  3 

[0./90.]  laminates.  The  difference  in  the  G  curves  for  these  two  laminates  is 
2  2  s 

greatest  at  small  values  of  crack  density.  Hence,  the  onset  of  90°  matrix  cracks 
may  occur  differently  for  the  two  laminates,  but  the  accumulation  of  cracks 
should  be  similar  for  both.  Figure  8  also  shows  a  comparison  of  (0/90) 
glass/epoxy  and  graphite/epoxy  laminates.  As  expected,  graphite/epoxy  shows  a 
lower  and  flatter  G  curve,  because  90°  ply  matrix  cracking  has  a  lesser  effect 
on  laminate  stiffness  in  graphite/epoxy  compared  to  glass/epoxy  [2J. 


Figure  9  compares  strain  energy  release  rate  as  a  function  of  crack 
density,  for  a  [0/90]  glass/epoxy  laminate,  calculated  both  from  eq.  8,  and 
from  eq.  7  using  dE/dD  derived  from  finite  element  results.  The  difference  in 
the  G  values  from  the  two  techniques  results  from  the  difference  in  slope  of  the 
stiffness  vs.  crack  density  expressions,  which  can  be  observed  in  figure  5.  At 
the  moment,  it  is  not  possible  to  establish  which  of  the  two  models  i3  closer  to 
the  actual  behavior  in  the  initial  region  of  low  crack  density,  since  both  rely 
on  some  simplifying  assumptions,  and  the  scatter  in  the  experimental  data  shown 
in  fig.  5  makes  it  difficult  to  resolve  which  analysis  i3  more  appropriate.  In 
the  high  crack  density  region,  however,  the  finite  element  analysis  seem3  to 
comply  better  with  the  experimental  data  and  with  the  assumption  that  for  high 
crack  density  G  must  reach  zero  in  order  to  explain  the  observed  phenomenon  of 
crack  density  saturation. 

Matrix  crack  onset 


To  evaluate  G  for  the  onset  of  90°  ply  matrix  cracking  in  a  [0  /90  ] 
c  m  n  s 

laminate,  it  is  easier  to  use  a  simpler  form  of  eq.  8  with  D  set  to  zero 


c  m  Ej  X  c 


A  similar  expression  was  derived  from  a  shear  lag  analysis  in  ref. [10]  and  [11]. 
Equation  9  differs  only  in  that  the  factor  3  included  in  the  X  expression  (eq.1) 
is  replaced  by  a  4  in  ref. [10]  and  [11].  This  difference  is  a  result  of  the 
different  distributions  of  the  longitudinal  displacements  in  the  transverse  ply 


used  in  the  shear  lag  analysis.  This  distribution  was  parabolic  in  the 
derivation  of  eq.  1  [13]  and  linear  in  ref.  10. 


For  a  generic  laminate  with  90°  plies  surrounded  by  0°  plies,  dE/dD  is 
easily  calculated  from  eq.(2)  when  D=0.  This  calculation  introduces  an 
approximate  linearized  expression  for  effective  damaged  lamina  transverse 
modulus,  obtained  from  eq.(1),  as  D  approaches  zero. 

Crack  growth 

As  crack  density  increases,  eq.  8  gives  a  lower  value  of  G  for  a  fixed 
strain.  In  a  constant  strain  fatigue  test  no  more  cracks  should  appear,  as  soon 
as  G  is  reduced  below  its  critical  value.  In  the  present  study,  where  all  the 
tests  were  run  under  load  control,  an  experimental  proof  was  not  possible  and, 
as  already  stated,  in  most  cases  the  matrix  cracks  reached  a  saturation  density 
value.  The  presence  of  this  saturation  density  is  not  totally  explained  with  the 
closed  form  expression,  because  in  eq.8,  the  dE/dA  term  never  goes  to  zero.  If 
dE/dA  does  not  go  to  zero,  then  an  increased  stress  will  always  generate  more 
cracks.  The  finite  elements  analysis  however,  shows  more  clearly  that  dE/dA  is 
very  close  to  zero  when  the  crack  density  is  near  to  the  saturation  value,  no 
matter  what  the  stress  level  (fig.  9). 

Looking  at  the  problem  from  another  point  of  view,  the  "apparent"  critical 
G  value,  GD ,  would  increase  with  increasing  crack  density.  Therefore,  a  crack 
resistance  curve  concept  was  borrowed  from  classical  fracture  mechanics  [17], 
showing  the  progressive  laminate  resistance  to  formation  of  new  matrix  cracks  as 
the  number  of  these  cracks  increases.  This  increased  crack  growth  resistance  is 
not  attributed  to  other  energy  absorbing  damage  mechanisms,  but  results  only 
from  the  constraint  and  load  transfer  capability  of  the  neighboring  plies 
adjacent  to  the  cracked  ply. 


Crack  growth  resistance  curves  (R-curves)  for  90°  ply  matrix  cracking  were 
obtained  from  quasi-static  tension  tests.  Five  tests  of  this  kind  were  performed 
on  [0/90]  and  [±45/0/90]  laminates.  Figure  10  shows  Gn  calculated  from  strain 
and  crack  density  data  for  both  laminates,  using  equation  (8).  The  first  point 
plotted  for  each  curve  corresponds  to  the  onset  of  the  first  matrix  crack,  and 
shows  Gc  as  would  be  calculated  from  eq.  9.  The  large  scatter  observed  in  quasi¬ 
isotropic  laminate  data  for  high  values  of  crack  density  is  mainly  due  to  the 
superposition  of  other  damage  modes  (i.e.  delamination)  that  can  influence  the 
laminate  behavior.  Results  from  the  two  laminates  are  superimposed  in  fig.  10 
and  show  good  agreement,  thus  indicating  that  the  R-curve  approach  is  suitable 
for  predicting  matrix  crack  growth  in  90°  plies  of  a  laminate,  when  the  cracked 
90°  plies  are  surrounded  by  0°  plies  of  the  same  material. 

FREE  EDGE  DELAMINATION 

A  model  for  stiffness  loss  and  strain  energy  release  rate  for  free  edge 
delamination  onset  and  growth  in  graphite/epoxy  laminates  has  been  presented  by 
O'Brien  [2],  In  this  model,  stiffness  decreases  linearly  as  delamination  size 
increases  and  can  be  expressed  by 


where  EQ  is  the  laminate  modulus  before  delamination,  a/b  is  the  ratio  of 

* 

delamination  size  to  specimen  width,  and  E  is  the  modulus  of  the  laminate 
completely  delaminated  along  one  or  more  interfaces. 


* 

The  delaminated  modulus,  E  ,  can  be  calculated  as  the  sum  of  the  products 
of  the  moduli  and  thicknesses  of  the  sublaminate3  which  are  created  by  the 
delamination.  That  is 
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yM 

^i-1 


Ei  fci 


(11  ) 


where  t^  and  are  the  thickness  and  the  moduli  of  the  ith  3ublaminate  and  M  is 
the  total  number  of  sublaminates.  The  strain  energy  release  rate  for  edge 
delamination  can  be  evaluated  by  differentiating  eq.  10  with  respect  to  a  and 
substituting  into  eq.  4  as  shown  in  ref.  2, 


G  = 


* 


E  ) 


(12) 


The  critical  value  of  strain  energy  release  rate  Gc  for  edge  delamination  is 
obtained  from  the  actual  onset  strain  .  This  approach  has  been  demonstrated  to 
be  useful  for  evaluating  and  predicting  free  edge  delamination  onset  and  growth 
in  graphite/epoxy  laminates  [14], 

For  glass/epoxy  laminates,  the  same  approach  has  been  attempted.  However, 
several  differences  were  noted.  In  graphite/epoxy,  stiffness  loss  due  to  90°  ply 
matrix  cracking  was  small  [2,6],  while  this  is  not  the  case  for  glass/epoxy 
laminates  [9-11,13].  In  addition,  because  G  for  90°  ply  matrix  cracking  is 
higher  for  glass/epoxy  than  graphite/epoxy  laminates  with  the  same  layup 

(fig. 8),  and  because  the  glass/epoxy  has  a  lower  Poisson  mismatch,  included  in 

* 

the  (Eq  -  E  )  term  of  eq.  12,  and  hence  a  lower  driving  force  for  edge 
delamination  than  graphite  epoxy,  the  glass/epoxy  laminate  is  more  likley  to 
have  extensive  matrix  cracking  before  the  onset  of  delamination.  These  matrix 


cracks  influence  the  stress  and  strain  field  in  the  region  where  the 


delamination  will  grow  and  they  act  as  a  delamination  resistance  mechanism  [16- 
20].  Also,  in  glass/epoxy,  because  of  the  lower  Poisson  ratio  mismatch,  edge 
delaminations  grow  only  a  little  in  quasi-static  conditions  before  the  final  0° 
fiber  failure.  Furthermore,  for  the  glass-epoxy  laminates  tested  in  fatigue,  the 
effect  of  delamination  on  laminate  stiffness  loss  is  partially  obscured  by  the 
simultaneous  growth  of  ±45°  matrix  cracking,  which  was  also  observed  previously 
for  graphite  epoxy  [20]. 

As  shown  earlier  in  fig.  3,  there  are  many  contributions  to  stiffness  loss 
in  a  [±45/0/90]s  glass/epoxy  specimen  during  a  fatigue  test.  To  evaluate  the 
contribution  of  each  damage  mode,  a  multiple  linear  regression  analysis  was 
performed  on  stiffness  loss  data.  Because  the  relationship  between  stiffness 
loss  and  edge  delamination  size  is  linear  (eq.10),  and  because  a  linear 
approximation  between  stiffness  loss  and  90°  ply  matrix  crack  density  seems 
reasonable  (fig. 5),  a  linear  relationship  was  assumed  between  stiffness  loss  and 
angle  ply  matrix  crack  density,  and  0/90  edge,  and  +45/-45  local  interface 
delamination  sizes,  as  follows: 

—  =  a0  -  aiD90  -  a2D45  -  a3D_45  -  a„(g)^  "  a^b^0/90  °3) 

The  set  of  coefficients  from  a0  to  a5,  for  which  equation  (13)  best 
approximated  the  experimental  data,  was  evaluated  statistically.  The  result  of 
this  analysis  is  shown  in  Table  1A.  As  previously  shown  in  fig. 3,  the  0/90 
interface  edge  delamination  and  the  angle  ply  cracking  occur  simultaneously. 
Therefore,  the  a2,  a3,  and  as  coefficients  are  given  by  the  regression  analysis 
with  a  high  degree  of  uncertainty.  The  contribution  of  each  damage  mode  to 
stiffness  loss  is  represented  by  the  product  of  the  regression  coefficient  in 
Table  1A  and  the  appropriate  damage  parameter  in  eq .  13-  Therefore,  the 


regression  parameters  alone  do  not  indicate  the  relative  contribution  of  each 
damage  mode  to  stiffness  loss.  For  example,  the  analysis  indicated  an  a2  value 
that  was  nearly  twice  a3,  but  fig. 3(a)  shows  that  the  crack  density  in  the  -45° 
plies  is  nearly  three  times  the  density  in  the  +45°  plies.  Hence,  the  external 
45°  matrix  cracks  had  less  influence  on  the  laminate  stiffness  loss  than  the 
internal  -45°  matrix  cracks.  Unfortunately,  there  is  currently  no  analytical 
model  that  directly  relates  stiffness  loss  to  angle  ply  matrix  crack  density  to 
verify  this  observation. 

In  order  to  achieve  a  better  experimental  evaluation  of  the  edge 
delamination  contribution  to  stiffness  loss  in  the  quasi-isotropic  laminate,  a 
similar  statistical  analysis  was  applied  to  the  static  test  data,  where  only  90° 
ply  cracks  and  edge  delaminations  were  found.  For  this  special  case,  eq.13 
reduces  to: 

f;  =  a0  -  alDgo  -  a8(|)  (14) 

The  set  of  coefficients  obtained  for  this  three  parameter  regression  analysis  is 
reported  in  Table  IB.  Dividing  eq.10  by  E0  and  comparing  to  eq.14  indicates  that 
a0  =  1 ,  and  that 

a5  =  (E0-  E*)/E0  (15) 

Therefore,  a  direct  comparison  may  be  made  between  a5  calculated  from  the 
regression  analysis  and  the  value  calculated  from  laminated  plate  theory  using 

eq.15.  According  to  the  regression  analysis,  the  delamination  contribution  to 

# 

the  longitudinal  modulus  loss  gives  an  (E0-E  )  value  of  <L 500  MPa.  However,  using 

* 

laminated  plate  theory  to  evaluate  the  terms  in  eq.15  yields  an  (E0-E  )  value 
ranging  from  10^0  MPa  to  770  MPa,  depending  whether  or  not  the  preexisting  90° 


p 


ply  matrix  cracks  are  considered.  Hence,  G  determined  using  (E0-E  )  from 
experiments,  where  significant  90°  ply  matrix  cracks  exist  ahead  of  the 
delamination  front,  would  be  larger  than  G  determined  from  calculated 
properties.  This  observation  is  consistent  with  previous  studies  that  showed  an 
edge  delamination  resistance  curve  could  be  generated  and  correlated  with  the 
increase  in  90°  ply  matrix  crack  density  L2,19,20j. 

These  results  suggest  that  the  interlaminar  stress  field  near  the 
delamination  front  in  the  0/90  interface  is  influenced  by  the  presence  of  90° 
ply  matrix  cracking.  A  three  dimensional  model  is  required  to  accurately  account 
for  the  interaction  between  matrix  cracks  and  delamination. 


FRACTURE  TOUGHNESS  EVALUATION 


From  the  analysis  above,  critical  strain  energy  release  rate  values  were 
obtained  for  both  matrix  cracking,  eq.(9),  and  edge  de lamina t  ion ,  eq.(12),  in 
static  tension  and  tension-tension  fatigue,  and  are  shown  in  Table  2.  In  both 
cases  dE/dA  values,  used  to  calculate  strain  energy  release  rate,  were  taken 
from  experimental  results. 

The  Gc  for  matrix  cracking  and  edge  delamination  both  decreased 

dramatically  in  fatigue  compared  to  the  static  values.  This  is  consistent  with 

similar  data  generated  on  graphite  reinforced  composites  with  a  variety  of 

matrices  [21-23J.  The  G  for  delamination  was  found  to  be  higher  than  the  G  for 

c  ^ 

matrix  cracking.  Although  different  fracture  modes  are  involved  for  these  two 
mechanisms,  this  mode  difference  is  not  sufficient  to  explain  the  large  observed 
difference  in  G  .  Other  sources  L 1 5  J  report  static  G  T  for  matrix  cracking  of 
reinforced  epoxy  composites  similar  to  the  measured  value  reported  in  Table  2. 
The  problem  of  evaluating  Gc  for  delamination  onset  is  again  closely  connected 


to  the  evaluation  of  the  stiffness  loss  due  to  delamination  and  is,  therefore, 
affected  by  the  uncertainties  reported  in  the  previous  section.  However  a 
comparison  with  data  previously  obtained  on  different  composites  materials  [16] 
showed  that  the  values  reported  in  Table  2  for  static  and  fatigue  onset  were  at 
least  congruent  with  those  data. 


! 
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CONCLUDING  REMARKS 


The  purpose  of  this  work  was  to  study  damage  onset,  accumulation  and  its 
consequences  on  glass/epoxy  multidirectional  laminates  using  simple  stress 
analysis  tools  and  fracture  mechanics.  The  relatively  lower  stiffness  of  glas3 
fibers  compared  to  graphite  fibers  and  the  higher  strains  involved  in  damage  and 
failure  processes  have  highlighted  the  considerable  importance  of  matrix  cracks 
parallel  to  the  fibers.  Glass-fiber  laminates  are  useful  in  studying  matrix 
cracking  but  they  show  incongruences  when  other  failure  modes,  such  as  edge 
delamination,  are  superimposed. 

A  simple  way  to  model  stiffness  reduction  due  to  cracked  90°  ply  matrix  was 
verified  and  shown  to  apply  over  a  wide  range  of  crack  densities  using  both 
F.E.M.  analysis  and  experiment. 

A  closed  form  expression  for  strain  energy  release  rate  associated  with 
matrix  cracking  has  been  proposed  and,  following  experiments,  strain  energy 
release  rate  resistance  curves  (R-curves)  for  two  different  layups  of  the  same 
material  have  been  shown  to  coincide  fairly  well.  The  strain  energy  release  rate 
approach  used  in  this  study  can  help  explain  the  observed  saturation  crack 
density  for  matrix  cracking. 

Edge  delamination  was  difficult  to  analyze  using  the  conventional  tools 
developed  for  graphite  reinforced  composites.  The  superpos i t ion  of  matrix 
cracking  leads  to  an  overestimation  of  the  interlaminar  fracture  toughness. 
Further  investigation  in  this  field  is  needed. 


APPENDIX  A  :  Finite  element  analysis  of  matrix  cracking 


Matrix  cracking  was  modeled  using  the  well  known  NASTRAN  finite  element 
computer  program,  both  to  obtain  a  comparison  of  eq.(1)  with  a  numerical  stress 
analysis,  and  to  evaluate  the  use  of  a  general  purpose  structural  program  in 
fracture  analysis  problems. 

The  shaded  area  of  fig.  4  was  modeled  for  different  crack  spacings.  The  only 
way  to  model  out  of  plane  fibers  with  the  NASTRAN  program  was  to  use  three- 
dimensional  elements.  A  single  thin  layer  was  modeled  using  8-noded 
isoparametr ic  elements  with  constraints  representing  a  plane-strain  situation. 


The  size  of  the  smallest  elements  near  the  crack  tip  was  1 / 1 8  of  the  ply 

thickness.  In  total  208  nodes  were  needed  for  the  [0/90]  laminate  model. 

s 

For  the  [ 0/90 J  laminate  the  prediction  of  laminate  plate  theory  was  exactly 
s 

matched  for  the  case  with  no  cracks  while  a  small  edge  effect  was  observed  in 

the  45°  and  -45°  plies  of  the  [±45/0/90]  laminate.  The  result  of  this  edge 

s 

effect  on  laminate  stiffness  was  eliminated  by  normalizing  each  stiffness  value 
obtained  for  different  crack  spacings  with  the  corresponding  value  obtained  for 
a  model  of  the  same  length  without  cracks. 


In  conclusion  it  has  been  demonstrated  that  a  general  purpose  program  like 
NASTRAN  is  suitable  to  model  this  problem,  although  the  use  of  three-dimensional 
elements  makes  it  rather  time  consuming  and  expensive. 
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APPENDIX  B  :  Exact  solution  for  longitudinal  modulus 


The  longitudinal  modulus  of  a  [0/90]  laminate  may  be  calculated  from 
classical  laminate  plate  theory  as 


E 


x 


(B1  ) 


and  due  to  the  symmetry  of  the  extensional  stiffness  matrix  [A]  of  a  [0/90J3 
laminate 


E 

x 


( B2) 


Introducing  the  reduced  stiffness  matrix  components  of  the  material  Q,lt  Q22, 
Q,  2  we  obtain 


E 

x 


\  [(Qi,  ♦  xQ22) 


Q,,2  ( 1 -x) 2 

CQx  x  +  Q22) 


(B3) 


where  x,  the  independent  variable  for  this  problem,  is  the  ratio  E2/(E2)0  of  the 
transverse  Young's  modulus  in  the  cracked  90°  ply  over  its  original  value. 

If  x=1  ,  eq.  B3  gives  E=E0,  the  stiffness  of  the  undamaged  laminate,  while  for 
x=0  it  gives  E=E^  used  in  eq.  3* 

Introducing  the  elastic  properties  of  T190/x7b1/b0  glass/epoxy,  the  maximum 
difference  between  eq.  B3  and  eq.  3  is  less  than  0.2  $  and  can  be  neglected. 
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Table  I.  Multiple  regression  parameters  for  fatigue  damaging  in  a  [±45/0/90j 

s 

glass/epoxy  laminate. 


A, 

.  FATIGUE 

TESTS 

—  =  a0  -  a!D90  -  a2D1J5 

a>%5  ' 

a*lg)  -a.(g) 

±45  0/90 

damage  type  parameter 

regression 

std.  dev 

undamaged  condition 

a0 

.9945 

.0042 

90°  ply  matrix  cracking 

ai 

.0674 

.0048 

+45°  ply  matrix  cracking 

a2 

.0415 

.0175 

-45°  ply  matrix  cracking 

a3 

.0238 

.0076 

+45/-45  interface  delamination 

a* 

.321 1 

.0302 

0/90  interface  delamination 

as 

.1488 

.0357 

Coefficient  of  determination 

=  0.89 

B. 

E 

E0  = 

.  STATIC 

a0  -  a,D 

TESTS 

90-  *.li) 

damage  type  parameter 

regression 

std. dev 

undamaged  condition 

a0 

1 .00 

.004 

90°  ply  matrix  cracking 

a, 

.11 

.005 

0/90  interface  delamination 

a5 

.10 

.039 

Coefficient  of  determination 

=  0.92 
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TABLE  II.  Measured  critical  strain  energy  release  rate  for  T190/x751/50 


damage 

90°  ply  matrix  cracking 
90°  ply  matrix  cracking 
0/90  interface  delamination 
0/90  interface  delamination 


load  condition 
static  tension 
tension-tension  fatigue 
static  tension 
tension-tension  fatigue 


Gc 

120  J/m2 
30  J/m2 
M60  J/m2 
110  J/m2 


^  "  •  t-  ■  n  ■  1 


EDGE  DELAMINATION 


MATRIX  CRACKS 
45°  PLY 


DEL AMINAT IONS 


LVDT  ROO 


SPECIMEN  EDGE 


Flg.1  :  Photograph  of  a  quasi  isotropic  L ±45/0/90 J  glass-epoxy  specimen,  during 
fatigue  test. 
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MATRIX  CRACK 


EDGE  DELAM 


CRACK  OENSITY  «  CRACK  WIDTH  :  D»w  [No#inm_1]  #  [nun] 


Stiffness  loss  in  static  and  fatigue  tests  on  [0  /90  ] 

tn  n 


glass/epoxy 


laminate. 


EQUATION  (1) 


EQUATION  (3) 


£  4  - 

o 
3 
I— 

HH  .3  - 

ID 


[0/90] g  LAMINATE 


[0/90] s  LAMINATE 


CRACK  0ENSITY*CRACK  WIDTH 
Fig. 6(a) 


90 3  PLY  EFFECTIVE  TRANSVERSE  MODULUS 
Fig. 6(b) 


CRACK  DENSITYaCRACK  WIDTH 
Fig. 6(c) 


Fig. 6  :  Evaluation  of  90°  ply  matrix  cracking  effect  on  generic  laminate 
stiffness. 


EQUATION  2) 


[I] 


CRACK  DENSITY  #  CRACK  WIDTH:  D*w  (No . amm-1] »  (mm] 


Strain  Energy  Release  Rate  versus  matrix  crack  density  for  different 
laminates . 


mmmm 


Comparison  of  Strain  Energy  Release  Rate  from  analytical  model  (eq.8) 
and  finite  element  analysis  in  a  [0/90]  glass-epoxy  laminate. 


▼  A 


A  A 


A  A 


▼  A 
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*\  ** 
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1+45/- 45/0/90 J g 


10/90] , 


i  I  I 
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Crack  density  *  crack  width.  Dw 


Fig. 10:  Matrix  cracking  resistance  curves  for  90°  ply  matrix  cracking  in  [0/90] 

and  [±45/0/90]  glass  epoxy  laminate, 
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